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Abstract 

The approach to the analysis of the dynamic of non-equihbrium open systems within the 
framework of the laws of classical mechanics on the example of a hard-disks is offered. This 
approach was based on Hamilton and Liouville generalized equations which was deduced for 
the subsystems of the nonequilibrium system. With the help of generalized Liouville equa- 
tion it was obtained that two types of dynamics are possible: reversible and irreversible. 
The connection between the dynamical parameter -generalized field of forces, and entropy 
is established. The estimation of characteristic time of establishment of equilibrium in the 
thermodynamic limit is realized. It is shown how from the condition of irreversibility of a 
hard-disk system, the condition of irreversibility for the rarefied system of potentially inter- 
acting particles follows. The explanation of the mechanism of irreversibility is submitted. 

Introduction. In connection with a modern physics, the Newton's equation is a background of 
dynamical picture of the world. It is caused by the fact that all four known fundamental interactions of 
elementary particles are potential. Potentiality of forces causes reversibility of the Newton equation in 
time. Hence, dynamics of all natural systems consisting of elementary particles should be reversible. But 
irreversibility is a basis of evolution. In fundamental physics irreversibility determines the contents of the 
second law of thermodynamics. According to this law there is function S named entropy, which can grow 
for isolated systems only, achieving a maximum in an equilibrium state. Thus, the fundamental physics 
includes areas contradicting each other: reversible classical mechanics and irreversible thermodynamics. 

So, the rigorous substantiation of the second law of thermodynamics within the framework of the 
classical mechanics laws is one of the primary tasks of modern physics. This problem keenly reveals diffi- 
culties and limits of capability resources of classical methods of physics. Despite of numerous attempts, 
this problem does not have satisfactory solution up to date. Thus, the question how one can describe 
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irreversible processes on the basis of the reversible equations of the classical mechanics, are discussed 
and studied [1, 5, 8, 13,14]. The solution of this question is necessary for development of the theory 
of irreversible processes on dynamic basis. In order to find an answer for this question, we studied the 
evolution of the non-equilibrium colliding hard-disk system [10-12]. 

Here we are going to show, how one can study of the mechanism of establishment of the equilibrium 
state within the framework of the classical mechanics laws if one replaces the canonical Hamilton and 
Liouville equations with the generalized equations, applied for nonholonomic and open systems. 

The main idea of our approach to analyzes of non-equilibrium system is splitting of conservative 
nonequilibrium system into interacting subsystems, research of dynamics of one of them with help of 
the generalized Liouville equation, is not brought to the assumption of potentiality of their forces of 
interaction. 

The work is constructed as follows. Using D'Alambcrt and the motion equations for disks, we obtain 
generalized Lagrange, Hamilton and Liouville equations for disk subsystem, which is selected from the 
whole disk system. Then the analysis of two types of evolution of nonequilibrium system - reversible 
and irreversible, was made. The existence of these types of evolution follows from the condition of 
compatibility of the generalized Liouville equation for subsystem and the canonical Liouville equation, 
applicable for whole system. Estimation of characteristic parameter of establishment of the equilibrium 
state was realized. Then, it is shown how the condition of irreversibility for the potentially interacting 
and rarefied disks system follows from the condition of irreversibility for a hard-disk system. The 
formula that determines the connection between increment entropy for full non-equilibrium system and 
the generalized field of forces of subsystems was obtained. In summary, the analysis of the results is 
submitted. 

Mathematical apparatus. Wc used an approach of pair interaction of disks. Their motion equations 
are deduced with help of the matrix of pair collisions. As complex plane this matrix is given [10] : 

(a —ib \ 
, where a = dkj exp{i'&kj)', f> = /3exp(ii?/jj); dkj = cos'&kj', P = sin'&kj', i- an 
—ib a I 

imaginary unit; k- and j- numbers of colliding disks; dkj- the impact parameter (IP), determined by 
distance between centers of colliding disks in the Cartesian plane system of coordinates with axes of x 
and y, in which the fc-disk swoops on the lying the j- disk along the x 

- axis. The scattering angle i^kj varies from to tt. In consequence of collision the transformation 
of disks velocities can be presented in such form: Vj^ — SkjVi~ (a), where Vj~ and V^"^ - are bivectors 
of velocities of k and j - disks before (— ), and after (-I-) collisions, correspondingly; Vkj={Vk,Vj}, 
Vj — Vjx + iVjy - are complex velocities of the incident disk and the disk - target with corresponding 
components to the x- and y- axes. The collisions are considered to be central, and friction is neglected. 
Masses and diameters of disks-" tf' are accepted to be equal to 1. Boundary conditions are given as either 
periodical or in form of hard walls. Prom (a) we can obtain equations for the change of velocities of 
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colliding disks [10-12]: 

Here, A^j =Vk — Vj- are relative velocities, dVk = — V^" , and SVj = Vj^ — - are changes of disks 
velocities in consequence of collisions, Lpkj = ifie^'^''^ ■ 

That is, Eq. (1) can be presented in the differential form of Somsikov, V.M. [2001]: 

Vk = ^kjS{i'kj{t))Akj (2) 

t 

where ipkj = [1 ~ Kfej|]/I^fejh S{tpkj)-delta function; lkj{t) = Zkj + J Akjdt - are distances between 



centers of colliding disks; z^j = ^k ~ ^j' ^'^'^ -^j ' ^'"^ initial values of disks coordinates; $fej = 
i{lkjAkj)/{\lkj\\Akj\). 

The Eq. (2) determines transformation of velocity of k -disk when it collides with j-disk. A right 
side of Eq.(2) represents force depending from Akj. Therefore the equation (2) is not a Newtonian one 
[4]. 

Let's take a disk system, which consist of N disks, when N ^ oo and — > oo, where N/LF' = finit; 
and Z/^ - is area occupied by disks. Then divide this system on R subsystems, so that in each subsystem 
will be T >> 1 disks. Therefore, N = RT. Let energy of the system is equal E = const. It is equal 
to the sum of internal energies of all subsystems and interaction energies between subsystems. Let 
us select one of them, which we call p-subsystem. Let us examine 6W^ - the virtual work of active 
forces made above p-subsystem. In the general case, this work can be presented us follows: SW^ = 

T N-T T 

J2 S ^ks^fk = J2 ^k^'f'k, where k = 1,2, 3...T -disks number of the p - subsystem, s = 1,2,3. ..N — T 

k=l s=l fe=l 

- external disk number to the p -subsystem of the disk, which interaction with k -disk of p-subsystem, 

N-T 

F^^ - the interaction force of fc - and s - disks, 5rk - virtual displacement of the k -disk, = Yl ^ks ■ 
Here we use that the virtual work of the interaction force of internal disks of a p-subsystem is equal to 
zero. 

In case of the pair interaction approaching, the virtual work of external forces gets the form: 5W^ = 

T T 

^ F^Svk = Yl ^ks^''"k, because in this case F^ = F^^. The inertial force can be presented so: SW-!^ = 

fe=l k=l 
T 

Y VkSrk. The sum of the active and the inertial forces is called the effective force. The principle of 
fe=i 

D'Alambert asserts that the work of effective forces is always equal to zero [2], i.e. 

6Wl = SW[^-5WS. (3) 
The feature above virtual work means, that in generally it is not reduced to a complete differential. 

R T 

Prom the motion equations for a hard disks it follows [11]: X] ( X] Pka)^'''^ = 0- Therefore the total 

p=i fe=i 

of the active, and the total inertial works for all subsystems at any moment of time is equal to zero, i.e. 

R R 

S ^^in ~ S — 0- This equaling can take place in two cases: when the sum of nonzero members 
is equal to zero, and when each member of the sum is equal to zero. It is obviously that the second case. 
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appropriate} to an equilibrium state, takes place when T — *• oo. For this case, with the help of the motion 
equations for a hard disks, it is possible to record: 

XT T 

E ^fc = E '^ksSiAs {t))Aks (t) = E ^ks = 0. (4) 

k=l fc=l k=l 

Equality to zero of the right-hand side of the equation (4) means, that the selected p-subsystem is 
in a stationary state. 

To obtain the general Lagrange equation for p-subsystem, let's transform D'Alambert equation (3) 
by multiplied it by dt, and integrated it over an interval from t = ti to t = t2- In the general case we 
have: JsW^dt^ J E - E ' = 

ti ti k—1 ji^k 

s'hi: V^dt - JiEiPP + E Fj^)Srk]dt - [E VkSrkYtl (b) In equation (b) the member |] F^^ 

ii k=l ti fe=l j^k k=l j^k 

determines the force of interaction a internal disks of the p-subsystcm. Then k and j -colliding disks 
from the p-subsystcm. The member Fj^ is the force on the p-s\ibsystem. If it is demanded that on the 
ends an interval [ti,t2] the virtual displacements are zero, the last member in (b) will be equal to zero. 

Let's neglected nonequilibrium inside a subsystem (this is a typical assumption for local equilib- 
rium). Then for internal forces of interaction of subsystem disks, we can set in the conformity such 
a function dependent on coordinates, ?7(ri, r2, ...ry), for which the following condition is satisfied: 

t2 T T t2 

J[J2 12 Flj5rk\dt = —S J U{ri,r2, ...rT)dt. Here - ri,r2...rT - coordinates of disks of the f^-subsystem 

ti k—l j^k ti 

disks. In general case it is impossible, to present acting on p-subsystem active forces, as a gradient of a 
force function [2]. In this case the equation (b) can be written as: 

= /E(i||- 1^21 -W* = (5) 

Ineqn. (5) we denote = E -2'+U{ri,r2, --rT)- Therefore, if the interaction of disks will be potential, 
the Lp wiU include also internal potential energy of the p-subsystem - U{ri,r2, ■■■Tt)- Because for any 
variations integral in equation (5) will be equal to zero, the next expression is carried out: 

T 

We denote T, K = ^p- 

k=l 

The equation (6) is a generalized equation of Lagrange for a p-subsystcm. So, Fp, is the polygenic 
force acting on the p-subsystem which dependent from its dynamic. When the Fp = 0, the equation (6) 
transforms to a canonical equation of Lagrange for an equilibrium, conservative system. 



Let us derive Hamilton's equations for p-subsystcm. The differential for Lp can be written as: 
.^dn + ^dVk) + ^dt, where |^ 



T 

dLp = Yl ("a^i^^fc + ffy'dVk) + -ofdt, where = pfe - is disks momentum. With the help of 



T T 

Lagrange transformation, it is possible to get: d[J2 PkVk — Lp] = '}2{—-g:f-drk + Vkdpk) gi''^^- 

k=i fe=i 



T 

Because -gf- = — -gf-, where Hp = PkVk — Lp, we will have from (6): 

fe=i 

f^'-A + K (7) 

These are the general Hamilton equations for the selected p-subsystem. The external forces, which acted 
on p-subsystem, presented in a right-hand side an equation (7). 

Using equations (7,8), we can find the Liouville equation for p-subsystem. For this purpose, let's to 
take a generalized current vector - Jp = {rk,Pk) of the p-subsystem in a phase space [13]. From equations 
(7,8), we find: 



The differential form of the particle number conservation law in the subsystem is a continuity equa- 
tion: ^ + div{Jpfp) = 0, where fp = fp{rk,Pk,t)- the normalized distribution function of disks in the 
p-subsystem. With the help of the continuity equation and equation (9) for a divergence of a generalized 

T 

current vector in a phase space, we can get: = i^k^+Pk§^) = %-+div{Jpfp)-fpdivJp = 

fe=i 

T 

-fp E ^Fk- So, we have: 

k=l 



Equation (10) is a Liouville equation for p-subsystem. It has a formal solution: 

d 
dpk ' 



= const -exph/d: ^Fl)dt\. 



fc=l 

The equation (10) is obtained from the common reasons. It is suitable for any interaction forces 
of subsystems. Thus, the equation (10) is applicable to analyze any nonequilibrium open systems. In 

T 

particularly, it can be used for explanation of irreversibility. The right side of (10), ^ -^F^, is the 
integral of collisions. This integral can be obtained from the motion equations of the systems element. 
For example, for a hard disks system it can be found with the help of the Eq. (2). 

The stationary nonequilibrium state of the system can be supported by external constant stream of 

T 

energy. Because ^ = ^ + J2 i^k^r- + Pk§^), then from (10) for stationary case we have: 

fe=i 

The equation (11) is similarly to the equation, which described fiuctuations of the distribution 
function in nonequilibrium gas [6]: 



Vk^=Stfp (12) 



Here, Stfp- is integral of collisions. 
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Let's consider the important interrelation between descriptions of dynamics of separate subsystems 

R T 

and dynamics of system as a whole. As the expression, J2 S -^fc = 0) is carried out, the next equation 

p=i fe=i 

for the full system Lagrangian, La, will have a place: 

dtdV, dr^-^ 

and the appropriate Liouvillc equation: + Vk^^ P^l^ ^ ^ T\ic fmiction, /tj, corresponds to the 

B. 

full system. The full system is conservative. Therefore, we have: ^ divJp = 0. This expression is 

p=i 

R R R — ^ R R 

equivalent to the next equality: |(E In/p) = |(ln U fp) = iU fp) TtiU fv) = 0- So, U fp = 

p=i p—i p—i p—i p—i 

R R 

const. In an equilibrium state we have Y\ fp = fn- Because the equality ^ i^p = is fulfilled during 

p=i p=i 

R 

all time, we have that equality, Y\ fp = fa, is & motion integral. It is in agreement with Liouville 
theorem about conservation of phase space [4] . So, only in two cases the Liouville equation for the whole 
system is in agreement with the general Liouville equation for selected subsystems: if the condition 

t T ^ T 

I m~^k)'^i ^ const (c) is satisfied when f ^ oo, or when, 'Sr^k)^ '^^ ^ periodic function 

k=l fe=l 

of time. The first case corresponds to the irreversible dynamics, and the second case corresponds to 
reversible dynamics. So, non-potentiality of interaction forces of subsystems causes in an opportunity of 
existence of irreversible dynamics. Let's consider, how these two types of dynamics can be realized. 

Irreversible dynamics. The condition (c) occurs if the next condition is satisfied: ^-Fp < 0, i.e., 
Fp{0,ojo) > Fp{t,uj) = S*Fp{0,ujf)). These conditions determine irreversibility. Here 5* - a operator of 
evolution or a phase stream; w - a point of phase space. Let's show that performance of these conditions 
follows from the law of momentum conservation of colliding disks. Firstly, we consider a simple case. Let 
us have P pairs of colliding disks. At that, P of them are resting disks, and P disks swoop on the latter 
with equal velocities directed along axis x. Let distances between colliding disks are equal, therefore all 
collisions occur through equal interval of time, r. In this case the total relative velocity of disks along x 

T 

axis before collision would be equal to PAk{t). After collisions the inequality PAk{t) > J2 ^kit+r) will 

fe=i 

take place, where Ak{t) are projections of relative velocity of disks on x axis. Equality would take place 
only if all collisions were frontal. So, as it follows from equation (2) the force may only decrease because 

T 

Fp is proportional to ^ Afe. Decreasing of Fp would take place in general case because inequality 
fe=i 

Fp ^0 takes place when the subsystem has average velocity which distinct from the average velocities of 

T 

external disks. As a result, the system would come to stationary state and the condition for it, ^ A^, 

fe=i 

would be satisfied. It is corresponds to the stationary state when Fp = 0. Prom here, the conclusion 
that the decrease of force, Fp, comes to homogeneous distribution of system energy, is follows. 

Let's consider the stability of the equilibrium state on the example of hard-disks system, and estimate 
characteristic time of decreasing force, -Fp. For this purpose we analyze solution's stability Eq.2 in the 
equilibrium point. 

The evolution of the p-subsystem is determined by the vector-column, V^, which components is 
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speed of disks of the p-subsystem: = {V^}, k = 1,2, 3...T. Some of the evolution's properties of this 
subsystem will be determined by the studying of the sum of it components. Let us designate this sum 
as Tp. Carrying out the summation in Eq.(2) on all disks of the p-subsystem, we shall obtain: 

T 

tp = J2^ksAks = Fp (14) 

The equation (14) describes the change of a total momentum, acting onto the p-subsystem as a result of 
collisions. The relaxation of a total momentum to zero is equivalent to relaxation to zero of the force, Fp. 
Now let us show, if the mixing property for a disks system is carried out, the homogeneous distribution 
of impact parameters of disks have place as well. In accordance to definition of the mixing condition, we 
have [7]: /i(5)//x(d) = 6/d where, n{d), is a measure corresponding to the total value of impact parameter 
- "(f'; 5 - is an arbitrary interval of the impact parameter and, ^{S), is a corresponding measure. The 
fulfillment of the mixing condition means the proportionality between the number of collisions of disks, 
falling at the interval, "S", and the length of this interval. It allows to say that distribution of the 
impact parameters is homogeneous. As it is well known [7, 13], for mixed systems the condition of 
depletion of correlations have place. For the equation (14) this condition can be written down so: 
< ^ks^ks >=< ^ks >< ^ks > i-e. average from two multiplied functions is equal to multiplication 
of these functions average. The ^ks is dependent from impact parameters, and Aks is dependent from 
relative velocities of colliding disks. Therefore this condition is similar to a condition of independence of 
coordinates and momenta that is widely used in the statistical physics Landau [1976] . Thus, it is possible 
to execute summation in the multiplier, ^ks, on impact parameters, independent from summation of 
expression 

Aks on relative speeds of colliding disks. Then under the condition of the homogeneous distribution 
of impact parameters and when T » oo, we can transit from summation to integration. So, we will 

T 1 

have [10,11]: (p = 1/T lim Yl ^ks = h I <Pksd{cos'&) = — |, where G = 2 is the normalization factor. 

k=l -1 

Taking it into account, we will have from eq. (14): 




The negative factor in the right side equation (15) means, that the force, Fp, will decrease. The stability 
of a stationary point p-subsystem can be established with the help of the motion equation (2). Let the 
point, Zq, be a stationary point, in which the, Fp, effecting on p-subsystem, is equal to zero. Prom the 
Lyapunov's theorem about stability follows that the point, Zq, is asymptotically stable if any deviation 
from it will be attenuated. Let us expand the left and right sides of the equation (16) into series by small 
parameter, v, of perturbation of velocities of disks of the p-subsystem, near point, Zq, and keep terms to 

T 

first order. The expansion of the left side of the equation (16) gives: i) = J2 ^k, where the summation 

k=l 

is carried out on components of the variation, v. In the expansion of the right side, the only remaining 

T 

part is, — |^efe = — 2/3t;. A contribution into the expansion will be given by collisions of disks of 

fc=i 
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the p-subsystem, with disks of its complement. So, we have: v = — |v- This equation means, that any 
system deviation from an equihbrium state will decries. Hence, the stationary point at performance of 
a mixing condition is steady. Stability is provided by occurrence of returning force. Fp, at a deviation 
of a subsystem from an equilibrium point. We shall note, that the equation (16) also follows from the 
theory of fluctuations (see Eq. [12]) where it is proved from other facts [3]. 

Reversibility dynamics. According to (c) and to the general concept of reversible dynamics the 
following inequalities should take place [7]: Fp{t + nTo,w) = Fp(t, w), AFq = S'*(t + nTo)Aro(t), where 
AFo is an element of volume of the phase space occupied with the subsystem at the moment, t; tq 
is period of system's return in the initial point; n = 1,2,3...; uj € AFq. This condition takes place 
when a vector of the force, Fp , rotates without change of its module. This is an example, demonstrating 
existence of such points, there is a system consisting of disks located on plane with a hard walls. If at the 
initial moment of time this system was at phase point, in which velocities of all disks are perpendicular 
to one of the wall, and all impacts of disks are frontal; so the condition of the reversibility would take 
place. I.e. in this case the system comes back periodically to the initial point. Therefore, such points 
are accepted to call as periodic or cyclic. Probability of system's return is determined by probability to 
be system in these points at the moment of preparation and by stability these points. Really, if at the 
initial moment of time the system is not in one of these points, the reversible dynamic is impossible and 
the system will never come to these points. Otherwise, it might leave this point, but that contradicts 
the definition of reversibility. 

It is possible to assume that for system of disks at TV — !■ oo, probability to appear in cyclic point is 
very small. Then, the probability of reversible dynamic is very small also, though it is vary from zero. 
Cyclic points are determined by symmetry of boundary conditions and groups of systems' symmetry, 
which can be found with help of matrix's collision. I.e. probability of reversibility of the system depends 
on geometrical characteristics of its elements and boimdary conditions. Research of these points, for 
example, their measures and stability, has basic interest for process of evolution [14]. So, we can say 
that reversibility is possible only at the event when at the moment of preparation of nonequilibrium 
system appears in one of cyclic points. Such interpretation of reversibility coincides with point of view 
of Einstein, according to which, probabilistic description of the system is dictated by probability to take 
one of the points in phase space in moment of the system preparation. But the dynamics of the system 
should be determined. 

I.e. the statistical description of systems is connected not with probabilistic the nature of processes, 
but with opportunity of use of such description for the analysis of dynamics of many-body systems. 

Irreversibility of rarefied system,s of potentially interacting disks. Let's consider, how , the irreversibil- 
ity for rarefied systems of potentially interacting disks is follow from the property of irreversibility for 
hard disks. Let us take into account, that pair collision approach is convenient for rarefied system. 
Dynamics of potentially interacting disks is described by the Newton's equation. In this equation the 
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force, Fks, between k and s disks is expressed by meaning of scalar function, U: Fks = where U is 

potential energy, rks is distance between disks. It is possible to take a characteristic radius of interaction 
disks, Rint for rarefied system, which is much less than length of free path, Ic, i.e. Ic ^ Rint- Disks can 
be considered to bo free when their distance up to the nearest neighbor disks is more than Rint- Let's 
consider scattering of two disks. 

In coordinate system of the centre of weights, character of scattering is determined by the formula of 
[4]: ipo = J^'"* , Here, ipo is a scattering angle; r is a distance between scattering disks; 



p is an impact parameter; rmin is a square root of the formula under radical; m is mass of disks; A^^ a 
velocity of swoop disk, which is equal to relative velocity of interaction disks in laboratory co-ordinate 
system. In coordinate system of the centre of weights m = 1/2 (the weight of a disk has been accepted 
to equal unit). 

It is easy to show that velocity of disks, after they abandon the area of interaction, is possible to 
be determined under formulas for hard-disks, if to make replacement in them at formula of velocity (see 
[4]): ^feg = Itt — 2ipo\/ 2, cos 6 ks = dks- Here, 6ks, dks, are an angles of scattering of the disk and an 
impact parameter accordingly. Thus, trajectories of disks under condition of r > Rint are possible to be 
determined with the help of formulas (2) without integration of the Newton equation of inner area of 
interaction. 

For rarefied system of disks the conditions, Ic/V tint, take place. Here, Unt = Rint/V is character- 
istic time during which the disks are in the interaction region, V is characteristic velocity of disks. In this 
case we can make transition, Rint 0. Such transition corresponds to approach of hard collisions. I.e. 
dynamics of the system is described by the equation (2) rigorously. It follows from here that dynamics 
of rarefied system can be studied both with the help of the equation of Newton and the equation (2). 
At Rint — > or tint — * the calculations results in both cases completely coincide. Hence, dynamics of 
rarefied system of potentially interacting particles may be described with help of motion equation (2). 
So, the rarefied system of potentially interacting particles, as well as system of hard colliding particles 
also can possess irreversible dynamics. 

Interrelation of the generalized field of force with entropy. If we know the generalized field of forces, 
it is possible to obtain deviation of entropy for given nonequilibrium condition of system from volume 
entropy of the equilibrium state. Really, let the nonequilibrium system come to equilibrium. Then work 
of the generalized field of forces would go on increase of internal energy of the system Landau, L.D. 
[1976]. In this case entropy deviation, A5, can be determined with the help of following expression 
Somsikov, V.M. [2003]: 



where Ep is kinetic energy of subsystem (all subsystems have T number of the disks). I.e. with help of 





(16) 
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equations (17) the dynamic parameter of system - generalized field of forces is connected with thermody- 
namic parameter -entropy. And as the system aspires to achieve a condition with the minimal internal 
work appropriate to the minimal value of a generalized field of forces, the equation (17) corresponds to 
the principle of the minimal entropy production. It is not difficult to obtain (for example when R = 2) 
that the entropy increment for all system may be only positive (though in some subsystems it can be 
negative). This conclusion corresponds to existing results connected with entropy change in the system 
at occurrence in them of some structures, for example, turbulent structure. So, Klimontovich, Yu.L. [1] 
with the help of the S-theorems, offered by him, proved decreasing of the entropy when formations of 
structures have a place. The basic difference of the formula (17) from the formula for entropy offered by 
Yu.L. Klimontovich, is that the deviation entropy, AS*, is determined by relative to the nonequilibrium 
condition, which accepted as "back-ground"; but in the formula (17) AS' is determined by deviation of 
entropy an nonequilibrium state from an equilibrium state. 

Conclusion. So, we shown an opportunity of explanation of irreversible dynamics on the basis of 
formalism of classical mechanics, if this formalism expands by inclusion to it of the generalized Hamilton 
and Liouville equations. Formally, irreversibility follows from the right side of the generalized Liouville 
equation. The nature of this equation is connected with appearing in the system of the generalized 
field of forces, when it deviates from equilibrium, and dependence its forces from particles velocities. 
This dependence is caused by that that moving particles create this field of forces. So one can submit 
the following explanation of irreversibility. Let the nonequilibrium system, which velocity of center of 
weights, Vo, is equal to zero. Because of nonequilibrium, velocity of the center of weights some subsystems 
will differ from zero. I.e. these subsystems will move relatively each other. Then the corresponding 
part of kinetic energy of the subsystems which connected with this relative subsystem moving will be 
redistributed between them proportionally to their relative velocity. This energy will go on increasing of 
system entropy. As a result the relative velocity of subsystems will go to zero and the system will come 
to an equilibrium state with maximal entropy. This is a physical essence of irreversibility. 

As in equilibrium state the generalized field of force is equal to zero, so in this case, the genciralized 
Hamilton and Liouville erjuations transform into the canonical equations. And therefore in equilibrium 
state Poincarc's recurrence theorem (sec. [13]) is applicable. Dynamic of such systems is completely 
reversible. But reversibility is possible not far from equilibrium state. So, it is because the factors, which 
determining irreversibility, has the second order of value in relation to linear disturbances. 

Thus, the analysis of evolution of open nonequilibrium systems becomes possible within the frame- 
work of the laws of classical mechanics only on the basis of the generalized Hamilton and Liouville 
equations. 
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